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High  Resolution  Stratospheric  Winds  From  Chemical 
Smoke  Trait  Experiments  at  White  Sands 
Missile  Range  and  Wallops  Island 

1.  INTRODUCTION 

During  the  past  several  years,  the  Air  Force  Geophysics  Laboratory  (AFGL). 
Aeronomy  Division  (LKD)  has  conducted  a  number  of  rocket  experiments  for  the 
purpose  of  releasing  in  a  controlled  manner  a  thick  trail  of  titanium  tetrachloride  to 
make  high  resolution  measurements  of  stratospheric  winds  and  windshears.  The 
measurements  show  jagged  wind  profiles  consistent  with  the  view  that  turbulence  in 
the  stratosphere  occurs  in  thin  layers  that  form  randomly.  There  are  structural 
features  in  the  wind  field,  corresponding  to  high  spatial  frequencies,  that  lead  to 
high  shears  over  short  vertical  distances.  These,  in  turn,  are  responsible  for  the 
onset  of  instabilities  that  grow  into  flattened  turbulent  regions,  accounting  for  some 
of  the  vertical  transport  of  tracers  or  pollutants  injected  into  the  stratosphere. 

Zones  of  intermittent  mixing  are  thus  created  which,  in  the  course  of  time,  have  an 
effect  analogous  to  a  diffusion  process  with  a  characteristic  diffusion  parameter. 

The  present  report  is  concerned  with  a  general  description  of  improved  technique 
used  to  reduce  and  analyze  the  raw  photographic  data.  It  also  presents  the  wind  and 
windshear  data  gathered  in  the  course  of  five  rocket  flights  for  which  the  triangulation 
has  been  found  to  be  both  consistent  and  reproducible.  Finally,  the  report  contains 
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a  discussion  of  some  possible  means  of  improving  the  accuracy  and  resolution  of 
stratospheric  winds  computed  from  smoke  trail  data.  The  transport  problem  has 
been  addressed  by  Dewan.  1 

2.  COLLECTION  OF  RAW  DATA 


Chemical  vapor  or  smoke  trails  dispensed  from  rockets  have  been  used  for 
over  two  decades  to  serve  as  visible  tracers  of  atmospheric  wind  speed  and  direc¬ 
tion.  For  stratospheric  altitudes  AFGL/LKD  utilizes  a  mixture  of  titanium  tetra¬ 
chloride  and  water -methanol  solution  to  form  a  thick  trail  that  follows  the  motions 
of  the  wind  field,  and  remains  well  defined  for  periods  of  several  minutes.  At 
later  times  the  spreading  of  the  trail  may  be  used  to  estimate  the  magnitude  and 

altitude  dependence  of  the  molecular  diffusion  coefficient.  Details  of  payload  con- 

2  3 

figuration,  and  vapor  release  hardware  are  given  by  Vickery  and  by  Stokes  et  al. 
The  trails  are  photographed  at  regular  intervals  from  the  time  of  their  injection 
into  the  ambient  air  until  they  drift  out  of  the  field  of  view  of  the  cameras,  or  have 
become  too  tenuous  to  be  photographed.  In  most  cases  this  happens  some  fifteen 
minutes  after  the  chemical  smoke  is  vented  from  the  rocket.  Figure  1  shows  the 
evolution  of  one  trail  (Trail  "20  MAY")  photographed  from  a  triangulation  station 
located  about  25  km  northwest  of  the  launching  site.  The  cameras  are  located  at 
two  or  more  ground  stations  chosen,  when  possible,  to  optimize  the  geometrical 
factors  that  influence  the  accuracy  of  the  triangulation  procedure  that,  subsequently, 
will  reconstruct  the  successive  positions  of  each  point  on  the  trail.  We  shall  not 
dwell  here  on  the  characteristics  of  the  cameras,  the  film,  the  exposure  sequence 
or  the  film  development  needed  to  obtain  negatives  with  sufficiently  high  contrast 
and  low  granularity  to  allow  the  triangulation  to  have  a  vertical  resolution  of 

10  meters  or  better.  These  characteristics  were  listed  in  an  earlier  report  by 
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Quesada  and  Trowbridge.  In  the  following  sections  we  shall  describe  changes  and 
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improvements  introduced  into  the  measurement  of  the  trail  and  star  background 
images,  and  in  the  triangulation  procedure  since  publication  of  the  report  just 
referenced. 
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Figure  1.  Evolution  of  Trail  "20  MAY".  Launched  at  Wallops  Island  on 
20  May  1978,  and  photographed  from  Triangulation  site  near 
Pocomoke  City,  Maryland 


2.1  Camera  Orientation 

Precise  wind  measurements  depend  on  tne  accuracy  with  whicn  the  orientation 
of  tne  camera  is  determined.  To  observe  the  motion  of  segments  of  a  trail  at  a 
50-km  range  at  10-meter  intervals,  we  cannot  tolerate  angular  errors  that  exceed 
0.  005  degree.  Fortunately,  when  considerable  pointing  errors  exist,  their  pres¬ 
ence  becomes  apparent  during  tne  triangulation  and  velocity  determination.  Tney 
produce  abrupt  discontinuities  in  the  trail  position,  and  excessively  large  wind 
shears.  Pointing  errors  are  introduced  by  a  variety  of  causes;  for  example,  wind 
gusts  that  momentarily  jar  tne  cameras,  accidental  operator  contact  with  the  camera 
mount,  differential  expansion  of  tne  mount  due  to  uneven  nesting,  mis  identification 
of  dim  stars  in  the  calibration  frames,  film  stretching  due  to  malfunction  of  the 
film  processor,  and  so  on.  Optimization  procedures  may  be  applied  to  reduce 
errors  of  tais  type,  but  taey  generally  produce  sligntly  different  wind  fields,  and 


we  no  longer  use  them.  Vector  procedures  are  used  to  determine  the  camera 
orientation  and  focal  length  from  photographs  of  stars  in  the  field  of  view  of  the 
cameras.  They  are  described  in  various  reports.  '  The  vector  procedures  re¬ 
late  the  star  image  positions  in  the  film  plane  to  their  equatorial  coordinates 
(right  ascension  and  declination).  The  computations  require  knowing  the  time  at 
which  the  calibration  photograph  was  exposed  and  the  camera  site  location  on  the 
Earth's  surface  (latitude,  longitude,  and  height  above  mean  sea  level).  These  are 
taken  from  geodetic  survey  information  with  a  target  accuracy  of  10  meters.  Camera 
azimuth,  elevation  and  horizontal  tilt,  and  the  precise  focal  length  of  the  lens  are 
determined  with  standard  deviations  consistent  with  the  desired  triangulation 
accuracy.  For  the  cameras  employed  at  the  WSMR  and  Wallops  Island  experiments, 
standard  deviations  of  0.  005  degree  for  the  angles,  and  0.  01  cm  for  the  focal  lengths 
were  routinely  computed  when  12  to  15  stars  were  used  to  determine  the  camera 
parameters.  These  deviations,  translated  to  spatial  positions,  produce  errors  of 
the  order  of  5  meters  at  the  (typically)  50-km  range  that  separates  the  cameras 
from  the  center  section  of  a  stratospheric  smoke  trail. 

2.2  Film  Meafurements 

In  principle,  it  is  possible  to  use  pairs  of  photographic  negatives  taken  simul¬ 
taneously  from  two  separate  ground  stations  to  reconstruct  the  spatial  location  of 
the  trail  at  the  instant  it  was  photographed.  To  do  so,  we  must  know  the  geo¬ 
graphical  location  of  the  site  and  the  orientation  of  each  camera.  The  orientation 
of  the  cameras  is  usually  defined  by  giving  the  azimuths  and  elevations  of  the  optical 
axes.  In  addition,  we  must  know  the  focal  length  of  the  lenses.  In  practice  there 
are  many  pitfalls  which  may  introduce  small  but  unacceptable  errors  in  the  triangu¬ 
lation.  Some  have  been  discussed  by  Trowbridge,  1  together  with  means  to  minimize 
or  compensate  for  their  effects.  The  basic  principle  is  to  measure  very  accurately 
the  coordinates  of  closely  spaced  points  on  the  center-line  of  the  trail  image.  To 
establish  a  reference  coordinate  system  common  to  all  photographs,  fiducial  marks 
are  recorded  on  every  frame  when  the  trail  is  photographed.  Since  several  hundred 
to  a  few  thousand  coordinate  pairs  must  be  measured  on  each  frame,  a  semiautomatic, 
computer  controlled  densitometric  system  was  designed  and  built  to  measure, 
digitize,  and  record  these  coordinates  on  magnetic  tape.  The  instrument  and  the 
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software  needed  to  make  it  operatioi  ,1  were  described  by  Trowbridge  and  Andrus. 
Recently,  a  series  of  modifications  was  needed  to  reestablish  the  operational 
status  of  the  instrument,  many  of  whose  digital  components  were  over  10  years  old 
and  could  not  be  replaced  or  repaired  when  they  malfunctioned.  The  instrument 
is  a  video  based  system  with  an  encoded  stage  that  transports  the  film,  and  upon 
command  from  the  computer  aligns  the  optic  axis  of  the  video  camera  with  any 
point  of  the  photographic  negative  under  examination.  Electronics  for  processing 
the  input  signals  were  designed  to  provide  adjustable  operating  ranges  of  approxi¬ 
mately  0.  5,  1.  0,  1.  5,  and  3.  0  optical  density  units  (D).  These  ranges  may  be 
shifted  to  cover  any  portion  of  the  total  useful  range  (4D)  of  the  densitometer.  The 
system,  shown  diagrammatically  in  Figure  2,  is  unique  in  that  gross  alignment  of 
features  to  be  measured  may  be  accomplished  using  distances  and  densities  from 
the  off-axis  video.  Fine  positioning  and  measurement  of  the  densities  of  desired 
features  are  performed  on  the  video  axis  as  the  film  is  aligned  with  the  axis  by  the 
moving  stage.  Thus,  effects  of  field  curvature,  distortion,  and  system  photometric 
nonuniformity  (vignetting)  are  essentially  removed.  Final  coordinate  information 
is  obtained  only  from  the  extremely  accurate  stage  drive  and  encoders.  The  maxi¬ 
mum  error  is  less  than  10  fJtm  over  a  150-mm  total  displacement,  with  repeatability 
and  precision  better  than  2  pirn. 

The  digitization  procedure  begins  with  the  alignment  of  the  film  frame  on  the 
scanner,  using  two  fiducial  markers  to  set  one  axis  of  the  film  coordinate  system 
perpendicular  to  one  axis  of  the  scanning  system.  A  third  fiducial  is  then  used  to 
determine  one  point  of  the  second  axis  of  the  film  coordinate  net.  whose  center 
coincides  with  the  camera  optical  axis.  The  trail  digitization  proper  starts  with 
the  introduction  of  a  trail  "skeleton,  "  that  is,  a  series  of  guide  points  widely  spaced 
along  the  trail  image  which  are  entered  manually  using  a  joy-stick  control  to  trans¬ 
port  the  film  to  the  desired  positions.  These  points  provide  starting  and  ending 
points  for  the  curve  following  software  of  the  computer  system,  and  allow  both 
accurate  determination  of  the  coordinates  of  the  trail  axis,  and  motion  in  the  proper 
direction  along  the  trail,  with  minimal  operator  intervention.  Proper  direction 
along  the  trail  is  crucial,  particularly  at  late  times  when  very  frequently  the  trail 
image  shows  closed  loops  that  must  be  traversed  in  the  correct  direction,  in  order 
to  avoid  serious  discontinuities  in  the  triangulation. 

Star  calibration  photographs  require  long  exposures  (2  min)  and  each  star 
image  is  a  narrow  arc  segment  on  the  film,  resulting  from  the  Earth's  rotation  during 
the  exposure  time.  The  coordinates  of  the  center  point  of  each  track  (correspond¬ 
ing  to  a  time  midway  between  the  beginning  and  end  of  the  exposure)  are  found  by 


8.  Trowbridge,  C.A. ,  and  Andrus,  W.  S.  (1978)  An  Automated  Coordinate  Measuring 
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first  digitizing  and  storing  coordinates  and  densities  measured  on  the  axis  of  the 
arc.  Density  and  position  information  are  then  transferred  to  the  computer,  which 
identifies  the  end  points  of  the  axial  arc  and  calculates  the  position  of  its  center. 

Star  tracks  have  been  found  to  provide  star  locations  with  a  higher  precision  than 
obtained  from  the  point  images  produced  by  short  exposures. 

2.3  Trail  Position  and  Horizontal  Velocity  Measurements 

The  triangulation  program  uses  the  coordinates  of  the  trail  centerline  as  viewed 
from  two  sites,  the  site  geodetic  coordinates,  and  each  camera  orientation  which, 
as  explained  previously,  is  derived  from  the  star  calibration  photographs.  In 
principle,  this  is  sufficient  to  reconstruct  the  spatial  location  of  the  trail.  We  have 
implemented  the  vector-matrix  techniques  described  in  References  5  and  6  to  create 
a  software  package  that  improves  significantly  the  ease  and  speed  of  the  triangula- 

9 

tion.  In  practice,  computational  problems,  such  as  described  by  Trowbridge  et  al, 
may  degrade  the  accuracy  of  the  results  over  short  sections  of  the  trail  and  reduce 
the  vertical  resolution  by  a  factor  of  approximately  2.  One  reason  is  that  the  trails 
generally  have  no  features  that  can  be  uniquely  identified  on  photographs  taken 
simultaneously  from  widely  separated  locations.  Triangulation  is,  therefore,  per¬ 
formed  using  an  iterative  approach  that  minimizes  discrepancies  in  the  dihedral 
angle  computed  from  each  site  to  a  given  point  on  the  trail.  The  quantity  which  is 
minimized,  and  serves  as  a  figure  of  merit,  is  the  angular  mismatch  between  the 
normals  to  the  dihedral  planes,  each  of  which  is  defined  by  the  common  line  through 
the  two  observation  sites  and  the  lines-of-sight  to  the  points  that  we  wish  to  match. 

The  spatial  position  of  the  trail  (altitude,  latitude  and  longitude)  is  determined, 
assuming  a  spherical  Earth,  from  the  intersection  (or  near  intersection)  of  the 
line-of-sight  vectors  for  points  that  have  been  matched.  The  data  are  stored  for 
later  use  by  the  computer  program  that  calculates  the  horizontal  velocities  from 
trail  positions  corresponding  to  known  time  intervals. 

Improvements  in  the  triangulation  procedure  were  obtained  by  a  change  in  the 
dihedral  angle  tolerances  that  resulted  in  point  matching  with  a  much  higher  accuracy. 
With  tigher  tolerances,  about  50  percent  of  the  digitized  points  are  discarded  during 
the  triangulation.  This  means  that  the  altitude  resolution  finally  obtained  is  some¬ 
what  less  than  10  m,  although  the  raw  data  were  digitized  with  about  a  5-m  resolution. 
Measures  of  "goodness  of  fit"  for  the  triangulation  other  than  control  over  the  dihedral 
angle  mismatch,  such  as  average  film  plane  mismatch  and  average  separation  vector, 
have  also  been  improved  by  a  factor  of  10.  There  is  also  excellent  agreement  between 
the  velocity  profiles  calculated  from  triangulation  site  number  1  against  site 
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number  2,  or  the  reverse  pairing. 

9.  Trowbridge,  C.A.,  Kofsky,  I.  L. ,  and  Johnson,  R.  H,  (1978)  Recording  and 
Analysis  of  Optical  Data  From  Stratospheric  Dynamics  Experiments^ 


Further  study  of  the  change  in  the  triangulation  procedure,  that  is,  insisting 
that  the  dihedral  angular  mismatch  be  extremely  small,  has  shown  that  it  very 
effectively  preserves  trail  continuity.  The  larger  values  of  dihedral  angle  mis¬ 
match  previously  tolerated  permitted  an  accumulation  of  error  to  occur,  and  grow 
up  to  a  limit.  This  was  followed  by  several  consecutive  mismatches,  and  then 
again  by  error  buildup.  The  resulting  sawtooth  error  generated  small  discontinuities 
(and  occasionally  larger  ones)  that  proved  to  be  troublesome  when  the  shears  were 
computed.  Owing  to  the  smaller  percentage  of  matches  obtained  when  a  tight 
dihedral  angle  tolerance  is  introduced,  some  resolution  is  lost.  We  have  estimated 
that  in  order  to  maintain  an  average  10-m  vertical  resolution,  the  raw  data  must  be 
digitized  at  intervals  of  about  12  /am  (on  the  film  plane)  for  a  50-km  range.  The 
exact  value  depends  on  the  degree  of  distortion  of  the  trail  induced  by  the  wind  field. 
By  reducing  the  separation  between  consecutive  scans  to  12  /an,  we  will  substantially 
increase  the  time  required  to  digitize  each  frame.  A  possible  reduction  in  the 
digitization  workload  has  been  suggested  in  Reference  7.  It  requires  digitizing  the 
films  from  one  site  at  the  nominal  10-m  resolution,  and  the  other  site  at  a  higher 
resolution. 

With  introduction  of  the  change  in  the  triangulation  procedure,  the  major  source 
of  error  for  determining  trail  positions  is  the  precision  with  which  the  camera 
orientation  can  be  measured.  Better  measurements  can  be  made  if  one  records  the 
start  and  end  times  of  the  exposure  with  an  accuracy  of  0.  1  second.  An  increase 
in  the  number  of  stars  used  (currently  10  to  15)  and  a  finer  digitizing  interval 
(typically  20  |im)  could  also  contribute  to  more  precise  camera  orientation 
parameters. 

The  next  step  of  the  smoke  trail  method  is  to  determine  winds  from  the  point 
trail  positions.  The  trail  position  data  from  the  triangulation  are  interpolated  at 
equal  altitude  intervals  by  either  a  cubic  polynomial  or  a  cubic  spline.  Then,  a 
least-squares  analysis  of  position  versus  time  is  introduced  to  derive  average 
horizontal  velocities.  Ordinarily,  we  can  use  a  sequence  of  4  to  6  frames  spanning 
an  interval  of  about  2  to  3  minutes.  Minor  changes  were  made  to  the  velocity 
routines  to  correct  inconsistencies  on  the  assignment  of  the  time  separation  for 
sequential  trail  positions.  The  velocity  values  were  unaffected  by  the  change,  but 
the  uncertainty  associated  with  each  velocity  was  reduced  by  nearly  a  factor  of  2. 

3.  WIND  AND  WINDSHEAR  RESULTS 

The  wind  measurements  which  will  be  presented  in  this  section  were  made  from 
raw  data  collected  from  1973  to  1978.  In  Appendix  A,  we  list  the  trails  and  star 
calibrations  whose  images  were  digitized.  For  the  trails,  some  features  of  the 


smoke  release  are  indicated.  Some  ot  the  trails  consist  of  up  to  10  segments.  By 
breaking  up  the  trail  in  this  manner,  we  can  unambiguously  identify  common  points 
on  different  views  of  the  trail,  and  thus  facilitate  the  operation  of  the  triangulation 
package.  Although  the  rocket  vehicle  and  the  smoke  release  mechanism  after  1973 
were  selected  to  cover  the  altitude  range  15  to  50  km,  it  is  often  impossible  to  use 
the  entire  length  of  the  trail.  For  example,  sometimes  the  camera  viewing  angles 
are  optimum  for  the  lower  sections  of  the  trail,  but  are  much  less  adequate  for  the 
upper  portions,  whose  images  show  multiple  overlapping  of  segments  or  extreme 
foreshortening  that  obliterates  structural  details,  or  does  not  permit  triangula¬ 
tion  with  the  required  resolution. 

The  following  pages  contain  the  results  of  our  measurements  in  graphical 
form.  In  Appendix  B  we  list  a  typical  output  of  the  triangulation  and  wind  programs. 
A  complete  data  file  for  all  trails  is  available  on  magnetic  tape. 

For  the  winds,  we  have  computed  the  East-West,  North-South  components, 
and  the  resultant  wind  speed.  We  have  also  computed  some  hodographs.  In  the 
shear  diagrams  we  have  displayed  the  values  of 


as  a  function  of  altitude. 

We  note  that,  in  general,  a  hodograph  covering  an  altitude  interval  of  many 
kilometers  in  the  stratosphere,  is  usually  of  limited  value,  because  of  the  convoluted 
nature  of  the  curve  traced  by  the  tip  of  the  vector  wind.  This  is  well  illustrated  in 
Figure  20,  where  altitude  ranges  from  18  to  33  km.  On  the  other  hand,  if  the  alti¬ 
tude  interval  extends  over  a  few  kilometers,  as  in  Figure  19,  the  hodograph  can  be 
very  informative. 
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Figure  4.  Windshear  vs  Altitude  for  Trail  "FLORA1 
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Figure  7.  Wind  and  Windshear  for  Upper  Section  of  Trail  "IRIS" 


Figure  8.  Wind  and  Windshear  vs  Altitude  for  Lower  Section  of  Trail  "22  APRIL". 
Rocket  launched  from  WSMR  on  22  April  1977  at  1215  UT 


Figure  9.  Wind  and  Windshear  vs  Altitude  for  Midsection  of  Trail  "22  APRIL" 
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Figure  10.  Wind  and  Windshear  vs  Altitude  for  Upper  Section  of  Trail  "22  APRIL" 
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Figure  11.  Wind  and  Windshear  vs  Altitude  for  Lower  Section  of  Trail  "20  MAY". 
Rocket  launched  from  Wallops  Island  on  20  May  1978  at  0929  UT 
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Figure  21.  Hodograph  for  Trail  "20  MAY" 


4.  DISCUSSION 

Measurements  of  stratospheric  winds  and  windshears  with  high  vertical  resolu¬ 
tion  show  that  the  dynamic  structure  of  the  stratosphere  is  very  complex.  A  salient 
feature  is  the  presence,  in  nearly  every  case,  of  a  variety  of  motions  corresponding 
to  scales  of  tens  of  meters,  and  varying  in  relatively  short  times,  superimposed  on 
longer  time  and  spatial  motions  that  represent  the  response  to  such  factors  as  tidal 
forces  and  diurnal  heating.  The  fine  scale  motions  are  particularly  apparent  when 
a  comparison  is  made  of  simultaneous,  or  nearly  simultaneous,  soundings  of  the 
stratosphere  by  means  of  techniques  that  have  widely  different  vertical  resolution. 
Figures  22,  23.  and  24  show  smoke  trail  measurements  done  during  an  earlier  pro- 
grame  at  WSMR,  and  rawinsonde  profiles  at  essentially  the  same  time.  Lack  of 
total  agreement  between  the  rawinsonde  and  the  smoothed  smoke  trail  profiles  is  to 
be  expected  because  of  time  differentials,  and  the  impossibility  of  matching  the  paths 
of  the  rawinsonde  balloon  and  the  smoke  trail  rocket.  The  comparison  shows,  how¬ 
ever,  that  gross  structural  features  of  the  wind  field  persist  over  time  intervals  of 
minutes  and  horizontal  distances  of  kilometers. 
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Figure  24.  Superposition  of  Figures  22  and  23 
to  Compare  Altitude  Resolution 


One  dash  for  the  trail  of  20  May  1978  was  redigitized  using  a  raSt,#?  '.!ur.i 
of  12  4m  (rather  than  the  more  common  48  jum),  in  order  to  detefstfne  whether  the 
denser  positional  data  would  improve  the  accuracy  and  vertical  resolution  of  the 
winds.  When  the  triangulation  results  were  scrutinized,  it  was  clear  that  the  aver¬ 
age  altitude  increment  had  been  reduced  by  approximately  a  factor  of  3.  However, 
no  improvement  was  observed  in  the  standard  deviation  of  the  resultant  winds.  It 
was  evident  that  improper  pairing  of  points  from  each  site  was  responsible  for  dis¬ 
continuities  in  the  triangulation  that  led  to  spurious  shears.  Repeating  the  triangula¬ 
tion  with  a  tigher  tolerance  for  the  dihedral  angle  did  not  improve  the  results.  We 
have  concluded  that  the  appearance  of  artificially  high  shears  is  often  caused  by 
attempts  to  reduce  data  with  a  resolution  finer  than  the  natural  limit  imposed  by 
the  precision  with  which  the  camera  orientation  and  focal  length  are  computed.  For 
the  measurements  given  in  this  report,  this  limit  lies  between  10  and  15  meters. 

On  the  other  hand,  smoke  trails  afford  in  practice,  as  our  measurements  demonstrate, 
a  technique  for  the  determination  of  horizontal  winds  at  intervals  of  10  meters  with 
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a  standard  deviation  in  the  velocity  that  averages  50  cm  per  sec.  These  numbers 
apply  to  large  sections  of  the  trails.  There  are  sometimes  short  intervals  along 
the  trails  where  the  measurements  are  not  that  good.  Loss  of  vertical  resolution 
by  a  factor  of  2,  and  wind  velocities  with  errors  exceeding  1  m  per  sec  occur  when 
the  trail  is  embedded  in  a  wind  field  that  aligns  a  short  section  of  it  with  a  line -of  - 
sight  from  one  of  the  cameras.  In  such  cases,  there  is  a  one-to-many  point 
correspondence  for  the  images  recorded  at  each  of  two  sites,  and  vertical  resolu¬ 
tion  as  well  as  wind  velocity  are  correspondingly  degraded.  In  principle,  it  is 
possible  to  refine  the  measurements  of  the  trail  images  to  provide  sampling  of  the 
wind  field  at  vertical  intervals  of  2  meters,  and  have  the  granularity  of  the  film  set 
the  limit  for  the  fine  scale  structure.  To  realize  this  improvement,  however,  it 
would  be  necessary  to  increase  the  accuracy  of  the  camera  orientation  and  focal 
length  determination,  which  at  present,  limits  the  resolution  to  no  better  than  about 
10  meters.  Obtaining  higher  accuracy  for  the  camera  parameters  does  not  appear 
(insurmountable.  We  estimate,  for  example,  that  a  factor  of  2  in  vertical  resolution 
can  be  gained  by  increasing  the  accuracy  in  timing  the  beginning  and  end  of  a  star 
calibration  exposure  to  0.  1  sec,  and  reducing  the  error  in  synchronization  of 
cameras  at  different  sites  to  a  comparable  value.  At  present  both  of  these  timing 
errors  vary  between  1/2  and  1  sec. 

Further  analysis  of  the  data  presented  in  this  report  is  in  progress.  It  will 
contribute  to  the  study  of  statistical  properties  of  stratospheric  wind  fields,  and 
will  form  the  basis  for  forthcoming  reports. 
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Appendix  A 

of  Data  Digitized 


Summary  of  Data  Digitised 


Event  Name 


or  Date 

Frame 

Site 

FLORA 

14 

Seehorn 

4  Jun  73 

17 

Seehorn 

19 

Seehorn 

21 

Seehorn 

23 

Seehorn 

14 

Hotel 

17 

Hotel 

19 

Hotel 

21 

Hotel 

IRIS 

65 

Seehorn 

6  Jun  73 

66 

Seehorn 

67 

Seehorn 

68 

Seehorn 

Data 

Single  Trail 
13-19  km 


Single  Trail 
14-23  km 


Event  Name 
or  Date 


Frame 


Site 


Data 


22  Apr  77 


20  May  78 


20  May  78 


20  May  78 


22  May  78 


81 

T5 

82 

T5 

83 

T5 

84 

T5 

81 

Two  Buttes 

82 

Two  Buttes 

83 

Two  Buttes 

84 

Two  Buttes 

85 

Two  Buttes 

2 

Pocomoke 

3 

Pocomoke 

4 

Pocomoke 

5 

Pocomoke 

6 

Pocomoke 

7 

Pocomoke 

8 

Pocomoke 

9 

Pocomoke 

2 

Wachapreague 

3 

Wachapreague 

4 

Wachapreague 

5 

Wachapreague 

6 

Wachapreague 

7 

Wachapreague 

8 

Wachapreague 

9 

Wachapreague 

4 

Pocomoke 

5 

Pocomoke 

6 

JPocomoke 

7 

Pocomoke 

4 

Wachapreague 

5 

Wachapreague 

6 

Wachapreague 

7 

Wachapreague 

41 

Pocomoke 

42 

Pocomoke 

43 

Pocomoke 

44 

Pocomoke 

45 

Pocomoke 

46 

Pocomoke 

47 

Pocomoke 

41 

Wachapreague 

42 

Wachapreague 

43 

Wachapreague 

44 

Wachapreague 

45 

Wachapreague 

46 

Wachapreague 

47 

Wachapreague 

Dashed  Trail 
18-37  km 


Dashed  Trail 
16-44  km 


High  Resolution  Scans 
of  Dash  3  -  Maximum 
Increment  *  12  4m 


Dashed  Trail 
29-42  km 


Summary  of  Star  Frames  Digitized 


Event  Name 


or  Date 

Frame 

Site 

FLORA 

9 

Hotel 

4  Jun  73 

11 

See horn 

IRIS 

63 

Hotel 

6  Jun  73 

63 

See horn 

22  Apr  77 

99 

T5 

99 

Two  Buttes 

20  May  7  8 

1 

Wachapreague 

22  May  78 

20 

Wachapreague 

1 

Pocomoke 

Maximum  Error 
(Along  Star  Track) 

±  24  Mm 
±  20  Mm 
±  24  Mm 
±  20  Mm 


Appendix  B 

Output  of  Position  and  Velocity  Programs  for  Trail  “IRIS" 

The  following  pages  contain  a  fraction  of  the  output  of  the  position  program  and 
the  entire  output  of  the  velocity  program  for  trail  "IRIS.  " 

Part  I.  Sample  Output  of  the  Triangulation  Program 

The  listing  of  the  position  program  displays  the  following  information: 

A.  First  Entry  If  a  match  is  found  for  a  point  belonging  to  site  1,  it  is 

blank.  If  no  match  is  found  the  program  writes  a 
shortened  line  that  begins  with  the  words  "NO  MATCH" 
and  then  the  number  of  points  on  the  file  corresponding 
to  site  2  that  it  tried  for  a  match  "(#)  ATTEMPTS". 
This  is  followed  by  the  point  numbers  for  which  the 
best  match  could  be  established  when  the  program 
examined  in  succession  points  from  the  file  for  site  2. 
The  match,  however,  does  not  lead  to  a  dihedral  error 
within  tolerance,  and  is  rejected.  Next  the  program 
writes  the  separation  vector  (meters)  corresponding  to 
this  best  match  and  the  dihedral  angular  discrepancy 
in  radians  X  1,000,000. 

B.  Second  Entry  The  point  number  on  the  file  of  film  plane  coordinates 

corresponding  to  site  1  that  is  to  be  matched  to  some 
point  on  the  file  for  site  2. 


C.  Third  and  Fourth 
Entries 


D.  Fifth  Entry 


E.  Sixth,  Seventh 
and  Eighth 
Entries 

F.  Ninth  and  Tenth 
Entries 

G.  Eleventh  Entry 

H.  Twelfth  Entry 


I.  Thirteenth  Entry 


J.  Fourteenth  Entry 


K.  Fifteenth  Entry 


Site  numbers.  The  triangulation  program  uses  two 
sites  at  a  time,  but  if  more  than  two  stations  were 
active  and  provided  photographic  data  one  can  use,  for 
example,  site  1  and  site  2,  or,  site  1  and  site  3,  and 
so  on. 

Numerical  value  of  discrepancy  (centimeters)  when  the 
point  on  site  1  which  has  just  been  matched  is 
projected  on  the  film  plane  of  site  2. 

Latitude  and  longitude  (degrees)  and  altitude 
(kilometers)  above  mean  sea  level  for  the  point  on  the 
trail  for  which  a  match  has  just  been  established. 

Ranges  to  each  site  (kilometers)  for  the  point  on  the 
trail  whose  coordinates  are  given  in  E. 

Figure  of  merit  for  the  present  match.  It  represents 
the  dihedral  angular  error  (radians)  X  1.000,000. 
Number  of  points  for  file  2  that  were  tried  before 
a  match  is  found.  If  the  first  attempt  leads  to  a 
dihedral  angle  within  tolerance,  the  program  does  not 
jump  to  the  next  point,  but  rather  tries  a  few 
neighboring  points  from  site  2  to  determine  whether 
they  result  in  smaller  dihedral  discrepancies,  and 
then  takes  as  a  match  the  pairing  with  minimum 
dihedral  error. 

Number  of  the  point  on  file  from  site  2  that  matched 
the  point  on  file  from  site  1  whose  number  is  listed 
in  B. 

Magnitude  (meters)  of  closure  error  of  lines -of -sight 
from  the  2  stations,  that  is,  length  of  vector  perpendicu¬ 
lar  to  both  line -of -sight  vectors. 

Flag  used  to  identify  printing  sequence  selected  for 
each  line  of  output. 
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Part  II.  Sample  Output  of  Velocity  Program. 

The  columns  of  the  velocity  program  are  labeled,  with  the  exception  of 
columns  7  and  9,  that  list  the  errors  in  the  velocity  components  in  meters  per 
second  to  within  the  nearest  decimeter  per  second. 

The  difference  between  columns  2  and  1  represents  the  time  elapsed  (in  seconds) 
between  the  first  and  last  trail  position  files  used  to  compute  the  velocity. 

Column  3  is  the  altitude  in  kilometers. 

Columns  4  and  5  give  the  azimuth  and  magnitude  of  the  velocity  vector.  For 
example,  an  azimuth  of  90  degrees  and  magnitude  of  20  meters  would  correspond 
to  a  wind  of  20  m/sec  blowing  from  West  to  East. 

Columns  6  and  7  list  the  East -West  component  of  the  wind  and  its  standard 
error. 

Columns  8  and  9  give  the  North-South  component  of  the  wind  and  its  standard 
error. 

The  last  column  indicates  the  number  of  frames  that  were  used  to  compute  the 
wind,  that  is,  through  how  many  points  a  least-squares  straight  line  was  fitted, 
whose  slope  is  the  wind  at  that  altitude. 
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